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radiance due to reflectance. Twenty eight NOAA 9/10 satellite 
passes are analyzed. Twenty two of the passes are found to contain 
a total of 316 ship tracks which is significantly more than that 
expected by earlier ship track studies. An existing ship track 
detection algorithm is used to conduct a statistical comparison of 
ship track and non-ship track, or ambient pixel reflectance of the 
NOAA 9 and 10 AVHRR channels 1 (0.63 /xm) , 3 (3.7 jxm) , and 
4 (11 /xm) . The results of the statistical analysis confirm, as 
found in previous studies, that the ship track pixels displayed a 
significant increase in values for channels 1 and 3 and a very 
slight increase for channel 4. 



in 



TABLE OF CONTENTS 



I . INTRODUCTION 1 

A. BACKGROUND 1 

B . AEROSOLS 2 

C. CLOUD MICROPHYSICAL PROPERTIES AND RADIATION 3 

1. Absorption and Emission 4 

2. Scattering 4 

3 . Reflectance 5 

D. SHIP EXHAUST PRODUCED CLOUDS 6 

E. OBJECTIVES AND ORGANIZATION 8 

II. DATA PROCESSING AND ANALYSIS 9 

A. SATELLITE 9 

B. SENSOR 9 

C. SIGNAL PROCESSING AND CALIBRATION 10 

1. Channel l 11 

2 . Channel 4 11 

3 . Channel 3 12 

a. Channel 3 Radiance 12 

b. Channel 3 Reflectance 12 

D. SHIP TRACK GEOGRAPHICAL STUDY 13 

E. STATISTICAL STUDY 13 

III. RESULTS 26 

A. SHIP TRACK FREQUENCY AND GEOGRAPHICAL ANALYSIS 26 

B. STATISTICAL STUDY 30 

IV. CONCLUSIONS AND RECOMMENDATIONS 59 

LIST OF REFERENCES 61 



INITIAL DISTRIBUTION LIST 



63 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



LIST OF FIGURES 

Fig. 1. Channel 2 overview 14 

Fig. 2. Channel 3 near-infrared subscene 15 

Fig. 3. Channel l visible subscene 16 

Fig. 4. Channel 1 ship tracks 17 

Fig. 5. Channel 3 ship tracks 18 

Fig. 6. Channel 3 subscene with 8 by 8 overlay 21 

Fig. 7. Coakley/Morehead algorithm on Fig. 6 subscene 22 

Fig. 8. Fig. 7 subscene reprocessed 23 

Fig. 9. Scatter plot, Chnl reflectivity vs. Chn3 radiance.. 24 
Fig. 10. Scatter plot, Chnl reflectivity vs. Chn4 radiance.. 25 

Fig. 11. Channel 1 ship tracks, 1-19 July 1987 27 

Fig. 12. Channel 3 ship tracks, 1-19 July 1987 28 

Fig. 13. Chnl/Chn3 ship tracks combined, 1-19 July 1987. ..29 

Fig. 14. Scatter plot from Table 4 - File number 8 34 

Fig. 15. Scatter plot from Table 4 - File number 9 35 

Fig. 16. Scatter plot from Table 4 - File number 10 36 

Fig. 17. Scatter plot from Table 4 - File number 11 37 

Fig. 18. Scatter plot from Table 4 - File number 12 38 

Fig. 19. Scatter plot from Table 4 - File number 13 39 

Fig. 20. Scatter plot from Table 4 - File number 14 40 

Fig. 21. Scatter plot from Table 4 - File number 15 41 

Fig. 22. Scatter plot from Table 4 - File number 16 42 



v 



Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 



23. Scatter plot from Table 4 - File number 17 43 

24. Scatter plot from Table 4 - File number 18 44 

25. Scatter plot from Table 4 - File number 19 45 

26. Scatter plot from Table 4 - File number 20 46 

27. Scatter plot from Table 4 - File number 21 47 

28. Scatter plot from Table 4 - File number 22 48 

29. Scatter plot from Table 4 - File number 23 49 

30. Scatter plot from Table 4 - File number 24 50 

31. Scatter plot from Table 4 - File number 25 51 

32. Scatter plot from Table 4 - File number 26 52 

33. Scatter plot from Table 4 - File number 27 53 

34. Scatter plot from Table 4 - File number 28 54 

35. Scatter plot from Table 4 - File number 70 55 



vi 



ACKNOWLEDGEMENTS 



I would like to thank Dr. Philip A. Durkee for his untiring 
support and guidance during the extended length of this project, 
without which this thesis would not have been completed. 

I owe a debt of gratitude to the Naval Postgraduate School 
Meteorology Department for making the final phases of this report 
as uncomplicated as it was. 

I am especially grateful to my wife for her encouragement, her 
faith that I could complete the task and for her many prayers. 



• • 
Vll 



I. INTRODUCTION 



A. BACKGROUND 

Much attention has been given to anomalous cloud lines 
seen in visible and near-infrared satellite images. The 
source of these anomalous cloud lines were first identified by 
Conover (1969) as exhaust from ships transiting at sea and are 
often referred to as ship tracks. He showed that the most 
likely cause of ship track clouds is the increase in the 
aerosol content produced by ship exhaust. Ship exhaust is a 
source of aerosols that act as cloud condensation nuclei (CCN) 
that have the effect of increasing the number of cloud 
droplets while reducing the droplet size (Coakley et al., 
1987) . 

Advanced Very High Resolution Radiometer (AVHRR) data was 
used in the Coakley et al. (1987) study with Channel 1 
at 0.63 jim, Channel 3 at 3.7 jim, and Channel 4 at 11 jra* 
Coakley et al. suggested that the effects of aerosols on cloud 
reflectivity may have a much greater influence on the earth's 
albedo than that due to the direct scattering and absorption 
by aerosols alone. They developed an algorithm to detect ship 
tracks automatically by comparing reflectivities of 
contaminated clouds to those of similar non-contaminated, or 
ambient clouds. The Coakley et al. algorithm worked well in 
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areas of uniform cloud cover, but did not perform well in 
areas without ship tracks, cloud free areas or regions of 
transition from one cloud regime to another. 

Morehead (1988) evaluated the Coakley et al. algorithm 
incorporating improvements in its ability to detect ship 
tracks. The modified algorithm was used in this study to 
statistically analyze and compare the radiative properties of 
the ship tracks found in data gathered from 1 to 19 July 1987 
in the North Pacific Ocean Basin. Three hundred and sixteen 
ship tracks were observed in this study. The Morehead 
modified algorithm failed to identify all the ship tracks that 
were found by visual inspection. It also had a further draw 
back of identifying non-ship track segments as ship tracks. 
A manual filtering technique was developed for this thesis to 
ensure that the statistical analysis was conducted only on 
confirmed ship track segments. This process will be described 
in detail in Chapter II. An additional focus of this study 
was to identify and plot the geographical location of the 
observed ship tracks. This process will be explained in 
Chapter III. 

B. AER080L8 

Aerosols alter the atmosphere's ability to scatter and 
absorb solar radiation. They also affect the radiation budget 
through their effect on clouds. The ability of aerosols to 
effect clouds comes from their function as CCN which serve as 
anchors on which cloud droplets can form. Atmospheric 
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particle size distributions cover a large range. CCN have a 
radius of about 0.1 /xm, cloud droplets have an average radius 
of 10 jxm, and large cloud droplets have a radius of 100's 
of /xm. While CCN are too small on their own to have a 
significant interaction with incoming solar radiation, once 
cloud droplets form around them their interaction with 
incoming solar radiation is very much greater. 

Ships are a source of aerosols in the atmosphere. Exhaust 
expelled onto the atmosphere is in the form of water vapor, 
gaseous constituents and aerosols. Aerosols are of the size 
to act as CCN, as discussed above, and boost the CCN 
concentration far above the normal level. Twomey (1968, 1984) 
and Hindman et al. (1977) have described how the concentration 
and size distribution of CCN have a direct impact on the 
formation of cloud droplets. The new source of CCN from ship 
exhaust will cause a change in reflective properties of 
clouds. Twomey and Cocks (1982) show that an increase in 
droplet concentration causes a decrease in the size of the 
droplet. The shift in droplet size, along with an increase in 
concentration, will result in a significant increase in 
reflectance due to scattering. 

C. CLOUD MICROPHYSICAL PROPERTIES AMD RADIATION 

Electromagnetic energy reacts with the atmosphere either 
by scattering or absorption. Energy that is absorbed is 
converted to an increase in the temperature of the absorbing 
body. The atmosphere, having a temperature, can also emit 
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electromagnetic energy. These processes represent the 
interactions that take place in the atmosphere. The source of 
radiation reaching a satellite is then due to either 
reflection or emission. Reflection results from energy 
entering the atmosphere and undergoing single or multiple 
scattering events until it is directed back out the top of the 
atmosphere, primarily affecting short wave solar radiation. 
Energy which is absorbed is converted into an increase in 
temperature of the absorbing body which can subsequently 
undergo emission at a multitude of wavelengths including long 
wave frequencies. 

1. Absorption and Emission 

The atmosphere has a direct impact on satellite 
measured radiance due to absorption and emission. Clouds are 
composed primarily of liquid water, so absorptive properties 
of water are a primary concern. Liquid water absorbs 
electromagnetic energy at all wavelengths. The absorption is 
much greater at the infrared wavelengths than in the visible. 
This leads to greater emission at the longer infrared 
wavelengths. 

2 . Scattering 

The type of scattering interactions depend on the 
cross-sectional area of the particle with which the energy is 
interacting. When the scattering particle size and the 
wavelength of the incident energy are of the same order, the 
interaction is best described by the Mie scattering theory. 
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For the particles in this category, the scattering properties 
become less dependent on the wavelength and more dependent on 
particle size. Mie scattering has a greater percentage of 
scatter in the forward direction than to the side or backward 
direction. As the droplet radius increases, forward 
scattering increases and backward scattering decreases. 

3. Reflectance 

In this study we are considering reflectance, the 
measure of energy reflected, from cloud surfaces. At the 
wavelengths considered, the reflectance depends most heavily 
on Mie scattering processes. A useful measure of the degree 
of scattering can be found by examining the optical depth of 
the cloud layer. Twomey (1977) concluded that pollution (or 
ship stack exhaust in the case of this study) , by increasing 
cloud nucleus concentration, hence increasing the numbers of 
cloud droplets; leads to increased cloud optical depth which 
increases cloud reflectivity. Hunt (1972) described the phase 
scattering characteristics as well as the particle size 
effects on reflectance of cloud layers in the visible and 
infrared windows. He showed that the changes in cloud 
radiative properties at the visible and infrared wavelengths 
are consistent with changes in the size distribution of 
particles such that, as the radius of the cloud droplets 
decreases the reflectance increases. 
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D. 8HIP EXHAUST PRODUCED CLOUDS 



Ship tracks appear as curvilinear features at visible and 
near-infrared wavelengths in satellite imagery. They have a 
plume like nature, narrow at the source and spreading 
horizontally with distance from the source. Ship tracks can 
be several hundred kilometers in length and can last for days. 

It appears that there are two different types of ship 
tracks. The first are seen in the visible satellite imagery 
where anomalous cloud lines form in areas where cloud 
formation is suppressed. These may be produced as the ship 
enters an area where the atmosphere is experiencing a CCN 
deficit. As the CCN from the exhaust gases mix with the air, 
the deficit is reduced enough for clouds to form. Another 
theory is that the cloud free region is in an area where the 
cloud droplets are of sufficient size to precipitate out 
therefore keeping the area free of clouds. As the ship passes 
through the area, the addition of the exhaust gas CCN reduces 
the size of the droplets so that they are no longer large 
enough to precipitate out and thus the cloud forms 
(Albrecht, 1990) . 

The second type of ship track can be seen in the near- 
infrared (NIR) satellite imagery. These are generated in 
cloudy areas and are detectable in the NIR due to the shift to 
smaller size cloud droplet distribution from the addition of 
CCN. These are not seen in the visible because at 0.63 nm 
there is no absorption. Therefore, the reflectance becomes 
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primarily a function of the liquid water property and size 
distribution. At 3.7 jum, there is moderate absorption hence 
reflectance becomes primarily dependent on the particle 
radius. At 11 /xm, there is high absorption and the particle 
acts as a black body with no reflectance. 

Fett et al. (1979) indicated that the most common areas of 
ship track formation were in regions of closed cellular clouds 
with medium base heights. Morehead (1988) pointed out that 
other conditions which help to promote ship track formation 
include areas with saturated or supersaturated air near the 
top of the marine layer; areas where the air temperature is 
cooler than the sea-surface temperature; or areas where there 
is a minimum of vertical wind shear. Saturated air at the top 
of the boundary layer provides sufficient moisture to form 
cloud droplets in the presence of the increase of CCN provided 
by the ship exhaust. The warm sea surface temperature 
enhances track formation in the presence of a moist boundary 
layer by creating an environment which is more favorable to 
fog and haze. The effect of wind shear in the boundary layer 
would be to disperse the CCN provided by the exhaust resulting 
in poorly defined or nonexistent cloud lines. Ship track 
cloud formation will be further enhanced in areas of 
stratus /stratocumulus cloud topped boundary layers, eastern 
oceans, and high latitudes. 
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E. OBJECTIVES AND ORGANIZATION 



The goal of this thesis is to analyze ship track data over 
an extensive data base of AVHRR satellite imagery using the 
ship track detection algorithm developed by Coakley et al. 
(1987) and modified by Morehead (1988) . This study will have 
two basic objectives. One is to develop a summary of where 
geographically the ship tracks most often occurred during the 
study period. An additional focus is to conduct a statistical 
evaluation of confirmed ship tracks to compare the 
relationship and contrasts between AVHRR channel 1 
reflectivity and AVHRR channels 3 and 4 radiance for ship 
track areas versus non-ship track areas. 

Chapter II will discuss the data used in the study as well 
as the satellite and sensor used to obtain the data. Chapter 
II will discuss the analysis used in the two parts of this 
study, including the techniques used to manipulate the data 
for determining ship track geographical locations and the 
techniques used in the statistical study. The results of the 
analysis will be presented in Chapter III and Chapter IV will 
contain the conclusions which are drawn from the results and 
recommendations for future work. 
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II. DATA PROCESSING AND ANALYSIS 



The first part of Chapter II describes the satellite 
platform and sensor used to collect the data for this study. 
Then, the source and format of the data will be discussed as 
will the basic signal processing and calibration which is 
conducted prior to use by the ship track cloud retrieval 
algorithm. 

A. SATELLITE 

The platforms providing the data for this study were the 
NOAA 9 and 10 polar orbiting satellite series sponsored by the 
National Oceanic and Atmospheric Administration (NOAA) . The 
satellites fly sun- synchronous polar orbits at an altitude of 
about 525 nautical miles. The NOAA 9/10 series provided 
coverage for different times of the day. NOAA-9 satellite 
passes were generally around 2300 UTC (1500 local standard 
time) and are roughly parallel to the North American coast. 
NOAA- 10 passes were at 1600 UTC (0800 local standard time) and 
cut NE-SW across the central California coast. The data were 
collected by Scripps Satellite Oceanography facility in La 
Jolla, California. 

B. 8ENSOR 

Radiance data on the NOAA 9/10 satellites are collected by 
the Advanced Very High Resolution Radiometer (AVHRR) 
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instrument. This sensor provided a nadir resolution of 1km by 
lkm at the polar orbit altitude of 525 N mi. 

Table 1 shows the five wave bands recorded by the AVHRR 
sensor. Channel 1 is in the visible, channels 2 and 3 are in 



TABLE 1 

WAVEBANDS OF THE ADVANCED 
VERY HIGH RESOLUTION RADIOMETER (AVHRR) 


Channel 


Wavebands (pm) 


Center 

Frequency (Mm) 


1 


0.58 - 0.68 


0.63 


2 


0.73 - 1.10 


0.83 


3 


3.55 - 3.93 


3.7 


4 


10.30 - 11.30 


11.0 


5 


11.50 - 12.50 


12.0 



the near-infrared, and channels 4 and 5 are in the thermal 
infrared portion of the electromagnetic spectrum. 



Channel 2 was used to prepare the satellite pass overview 
that distinguished what areas to investigate for possible ship 
tracks. Channels 1, 3 and 4 were used in the ship track 
detection algorithm. 

C. SIGNAL PROCESSING AND CALIBRATION 

The development of cloud analysis algorithms as part of 
the International Satellite Cloud Climatology Project (ISCCP) 
motivated the collection of the data used in this study. All 
the data were part of the First ISCCP Regional Experiment 
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(FIRE) conducted in the summer of 1987. All the high density 
tapes used for the analysis had been previously navigated for 
the FIRE study. 

The original data came on high density tapes in the form 
of counts of radiant energy in each of the five channels 
listed above. The radiant energy data were calibrated and 
converted to geophysical parameters. The calibration of 
channels 1, 3 and 4 is discussed briefly below and is 
discussed in full by Allen (1987). 

1. Channel 1 

This is the visible channel, which was used in this 
study to search for areas where the potential for finding ship 
tracks was greatest. Since it is known that ship tracks are 
not likely to be found in cloud free regions in the visible, 
these areas could be avoided. The calibration of channel 1 
data is accomplished assuming a linear relationship between 
counts received by the sensor and the reflectance. The 
calibration is done in terms of albedo and the results are in 
units of percent of reflectance. 

2 . Channel 4 

Channel 4 imagery is used to help classify cloud 
masses. It is used in the Coakley algorithm to help locate 
non-ship track pixels with similar microphysical properties as 
nearby ship track pixels. Channel 4 data result mainly from 
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thermal emission. The data are converted to radiance 
measurements with units of W/m 2 -sr-cm' 1 using a linear 
correlation to counts of radiant energy. 

3. Channel 3 

The channel 3 near-infrared band data contains the most 
important information for ship track detection. Ship tracks 
show up very clearly in channel 3 while they may not be 
detected in other wavelengths. Channel 3 data is comprised of 
both reflectance and thermal emittance information for 
daylight passes. As discussed by Morehead (1988) , two types 
of channel 3 data are used: 

a. Channel 3 Radiance 

The first is a channel 3 radiance which contains 
both the reflectance and emission contributions, determined by 
a linear correlation between counts measured and radiance, 
with units of W/m 2 -sr-cm' 1 . 

b. Channel 3 Reflectance 

The second is a channel 3 reflectance signature 
estimating only that portion of the measured irradiance 
resulting from reflectance. 

The reflectance is derived from the channel 3 
radiance and the channel 4 radiance which represents the 
thermal emission. The channel 4 data are used to estimate 
emission in channel 3. This emission is subtracted from the 
total radiance leaving only a reflectance. 
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D. 8HIP TRACK GEOGRAPHICAL STUDY 



One focus of this analysis was to determine geographically 
where the ship tracks formed during the period of the study 
(1-19 July, 1987) . Each AVHRR data tape processed contained 
an overview image consisting of one satellite pass of the 
North Pacific Ocean basin (Figure 1) . Each pass was 
subjectively scanned, subscene by subscene (500km by 500km) , 
in both channel 3 and channel 1 to locate ship track clouds 
(Figures 2 and 3) . Once a ship track was found, latitude and 
longitude positions were traced along the ship track. The 
ship tracks were then plotted to show their location. 
Geographical diagrams were produced for each pass showing the 
position and track of all ship tracks found in both channel 1 
and channel 3. Figure 4 and Figure 5 show the geographical 
plots of ship tracks found in channel 1 and channel 3 
respectively, for the 7 July pass whose overview is shown in 
Figure 1. 

E. STATISTICAL STUDY 

The second part of this study involved a statistical 
comparison of channels 1, 3 and 4. An algorithm developed by 
Coakley et al. (1987) and modified by Morehead (1988) was 
applied to each subscene determined to contain ship tracks, as 
found in the geographical process above. The result of each 
subscene statistical analysis was compared to the respective 
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Figure 1. Channel 2 overview: 

NOAA-9 satellite pass AR6085 (2237 
UTC, 7 July 1987) . Displays the 
location of the subscene used as an 
example of the data processing of 
this thesis. 
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Figure 2. Channel 3 near*- infrared subscene: AR6085 (2237 UTC , 
7 July 1987) . 
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Figure 3. Channel l visible subscene: 

7 July 1987) . 



AR6085 (2237 UTC, 
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channel 3 image. A numbered 8 by 8 grid of 64 pixel squares 
was overlaid on both the statistical image and the channel 3 
image. The numbered squares (1 - 64) facilitated in selecting 
which pixel areas to filter out. Only the grid areas with 
confirmed ship track data were reprocessed, filtering out non- 
ship track statistical data. 

This process is shown for the NOAA-9 pass on 7 July 1987, 
in Figures 6, 7 and 8. Figure 6 shows a channel 3 subscene 
with ship tracks clearly visible as dark curvilinear lines 
with an 8 by 8 grid overlaid. Figure 7 shows the results of 
the Coakley/Morehead algorithm processed on the subscene. 
Figure 8 shows the reprocessed analysis with the non-ship 
track lines filtered out and the normal or ambient control 
points included. The filtering process was conducted to 
create a statistical data base representing as near as 
possible, only confirmed ship tracks. The final step in the 
analysis was to process grid points containing confirmed ship 
tracks through a comparison routine of ship track vs. ambient 
pixels, giving values of reflectance for channel 1 and 
radiance for channels 3 and 4. The results of this routine 
were presented in tabular and scatter plot format with 
comparisons among channels l, 3 and 4. Figures 9 and 10 are 
the scatter plots for the case presented in Figures 6 through 
8 above. Figure 9 shows the comparison of channel 1 
reflectivity to channel 3 radiance. Ambient pixels are 
represented by ( * ) s and enclosed by a solid ellipse. 
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Ship track pixels are represented by ( + )s and enclosed by a 
dashed ellipse. The increase in reflectivity in channel 1 and 
radiance for channel 3 for the ship track pixels can be 
readily seen in Figure 9. Figure 10 shows a comparison 
between channel 1 reflectivity and channel 4 radiance. It is 
noted that there is very little difference between the ambient 
and ship track pixel values in this comparison. This result 
will be discussed further in Chapter III. 
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Figure 6. Channel 3 subscene with 8 Z 

(2237 UTC, 7 July 1987). 



8 grid overlay: AR6085 
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Figure 7. Coakley/Morehead algorithm on Fig 6. subscene. 
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Figure 8. Fig. 7 subsoene reprocessed: After grids 

containing non-ship track analysis were 
filtered out. 
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SHIP TRACK SCATTER PLOT 




Figure 9. Scatter plot, Chn 1 reflectivity vs. Chn 3 

radiance: Ship track pixels (+) and ambient 

pixels ( • ) . 
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SHIP TRACK SCATTER PLOT 




11-yrxM RADIANCE 
(rnWim -2 sr -1 cm) 



Figure 10. Scatter plot Chn 1 reflectivity vs. Ctan 4 

radiance: Ship track pixels (+) and ambient 

pixels ( • ) . 
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m. RESULTS 



A. SHIP TRACK FREQUENCY AND GEOGRAPHICAL ANALYSIS 

Coakley et al. (1987) made an assessment that ship tracks 
would occur in only 5 to 10% of all orbital passes. During 
the course of this study, 28 NOAA passes (17 NOAA - 9, 
11 NOAA - 10) were analyzed. This encompassed 17 days from 1 
to 19 July 1987. During this period, 408 subscenes (500km by 
500km) were examined leading to the discovery of 316 ship 
tracks. Table 2 provides a summary of the ship tracks found 



TABLE 2 

RESULTS OF GEOGRAPHICAL ANALYSIS 



PASSES WITH SHIP TRACKS 



SATELLITE 


PASSES 


IN CHN 3 
ONLY 


IN CHNs 
1 & 3 


TOTAL 


NOAA- 9 


17 


3 


10 


13 (76%) 


NOAA- 10 


11 


5 


4 


9 (82%) 


TOTAL 

1 


28 


8 (29%) 


14 (50%) 


22 (79%) 



in the analysis. Only six of the passes did not contain ship 
tracks in either channel 1 or 3. This analysis found that 79% 
of all passes studied contained ship tracks, which is 
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dramatically higher than the Coakley et al. assessment. 
Figures 11, 12 and 13 display the plots of ship tracks found 
in channel 1, 3 and combined respectively for all the passes 
processed for this study. 
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B. STATISTICAL 8TUDY 



Ship track studies by Coakley et al. (1987) and Morehead 
(1988) were conducted on special cases with only a limited 
data set and provided a snapshot of the radiative properties 
of a few ship tracks. The analysis of this thesis was 
conducted over a much broader data base. As the geographical 
analysis was being conducted, all subscenes found to contain 
ship tracks were annotated for further analysis using the 
modified Coakley algorithm. Once the filtering process, 
discussed previously, was completed, a statistical data base 
was created that provided a data enriched view of the 
radiative properties of the ship tracks detected. 

One of the results from the algorithm was the degree of 
variability in the radiative and reflective properties of the 
ship tracks. Table 3 presents the maximum and minimum values 
of channel 1 reflectivity and channels 3 and 4 radiance for 



TABLE 3 




AMBIENT 


PIXELS 


SHIP TRACK PIXELS 




MAX 


MIN 


MAX 


MIN 


CHANNEL 1 


39.5942 


16.9879 


45.1522 


21.2278 


CHANNEL 3 


0.9252 


0.4462 


1.0548 


0.5400 


CHANNEL 4 


96.6457 


76.6413 


96.6786 


76.6425 | 
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the ambient and ship track pixels processed. Table 4 gives 
mean radiance values for the channel 1, 3 and 4 ambient and 
ship track pixels, for the 61 subscenes processed. When the 
maximum and minimum values shown in Table 3 are compared to 
the average of the means from Table 4 , one can see that for 
the ship track pixels, the spread from the average is greater 
than the spread from the average for the ambient pixels. This 
implies that the response of clouds to ship effects is not 
confined to one resulting reflectance value. 

Figures 14 through 35 represent some selected scatter 
plots from the data set. The diagrams plot distributions of 
channel 1 reflectivity as a function of channel 3 or channel 
4 radiance. Pixels contaminated by ship-stack exhaust are 
designated by (+) and those of randomly selected, nearby 
ambient fields of view, are designated by (•)• The solid 
(ambient control group) and dashed (ship track) ellipses 
represent fits at 2 standard deviations to the reflectivity 
and radiance distributions, which are taken to be joint 
Gaussian distributions. 

The study by Coakley et al. (1987) displayed a similar 
scatter plot depicting an overall increase in both channel 1 
reflectivity and channel 3 radiance for the selected subscene 
of that study. It can be seen in Figures 14 through 35 that 
the visible reflectance is not always increased for the ship 
track group. In fact, only 36 of the 61 files presented in 
Table 4, had increases in channel 1 reflectivity. Figures 15, 
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21 - 25, and 27 all show a decrease for the visible 
reflectance. The scatter plots graphically show the great 
deal of variability in both the ship track and control group 
values. For example, Figure 16 shows a range in channel 1 
ambient reflectance of 12 - 32% and in channel 3 ambient 
radiance of 0.5 - 0.75 mWm‘ 2 sr' , cm. Figure 16 shows a range in 
channel 1 ship track reflectance of 19 - 31% and in channel 3 
ship track reflectance of 0.71 - 0.94 mWm' 2 sr' , cm. Coakley 
et al., attributed this variability to the irregularity in the 
distribution and size of cloud droplets within the same cloud 
mass, which serves to obscure some of the increase in channel 
1 reflectivity. It is also shown in Figures 14 through 35, 
that the channel 3 near-infrared radiance values are higher in 
every case represented in this study, confirming similar 
observations made by Coakley et al. 

Table 4 shows the small variability of the channel 4 
radiance and also the almost negligible variation between ship 
track and control group radiance levels. Coakley et al., 
attributed this to the observation that at the channel 4 
wavelength, liquid water is a strong absorber, therefore the 
amount of scattering is negligibly small. Also, there are 
enough droplets so that marine stratocumulus emit like black 
bodies regardless of droplet size and concentrations. Thus, 
there is no change in thermal emission between contaminated 
and noncontaminated clouds. 
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SHIP TRACK SCATTER PLOT 
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Figure 14. 



Scatter plot from Table 4 - File 
number 8. AR6105 (2314 UTC, 13/7/87). 
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SHIP TRACK SCATTER PLOT 
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Figure 15. Scatter plot from Table 4 - File 

number 9. AR6105 (2314 UTC, 13/7/87). 
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SHIP TRACI- SCATTEF PLOT 
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Figure 16. Scatter plot from Table 4 - File 

number 10. AR6105 (2314 UTC, 13/7/87) . 
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SHIP TRACK SCATTER PLOT 




Figure 17. Scatter plot from Table 4 - File 

number 11. AR6105 (2314 UTC, 13/7/87). 
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SHIP TRACK SCATTER PLOT 
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Figure 18. Scatter plot from Table 4 - File 

number 12. AR6105 (2314 UTC, 13/7/87). 
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SHIF TRACK SCATTER PLOT 
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Figure 19. Scatter plot from Table 4 - File 

number 13. AR6105 (2314 UTC, 13/7/87). 
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Figure 20. Scatter plot from Table 4 - File 

number 14. AR6113 (1526 UTC, 15/7/87) . 
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SHIP TRACK SCATTER PLOT 
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Figure 21. 



Scatter plot from Table 4 - File 
number 15. AR6114 (2252 UTC, 15/7/87). 
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SHIP TRACI' SCATTER PLOT 




Figure 22. Scatter plot from Table 4 - File 

number 16. AR6114 (2252 UTC, 15/7/87). 
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SHIP TFuO SCATTER PLOT 




Figure 23. 



Scatter plot from Table 4 - File 
number 17. AR6114 (2252 UTC, 15/7/87). 
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SHIF TRACK SCATTER PLOT 
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Figure 24. Scatter plot from Table 4 - File 

number 18. AR6120 (2231 UTC, 17/7/87). 
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SHIP TRACK SCATTER PLOT 
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Figure 25. Scatter plot from Table 4 - File 

number 19. AR6120 (2231 UTC, 17/7/87). 
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SHIP TRACK SCATTER PLOT 
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Figure 26. Scatter plot from Table 4 - File 

number 20. AR6049 (1529 UTC, 1/7/87). 
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SHIF TRACK SCATTER PLOT 




Figure 27. Scatter plot from Table 4 - File 

number 21. AR6049 (1529 UTC, 1/7/87). 
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SHIP TRACK SCATTER PLOT 
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Figure 28. Scatter plot from Table 4 - File 

number 22. AR6049 (1529 UTC, 1/7/87). 
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SHIP TRACK SCATTEP PLOT 
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Figure 29. Scatter plot from Table 4 - File 

number 23. AR6049 (1529 UTC, 1/7/87). 
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SHIP TRACK SCATTER PLOT 




(mWm 2 sr ' cm) 



Figure 30. Scatter plot from Table 4 - File 

number 24. AR6049 (1529 UTC, 1/7/87). 
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Figure 31. Scatter plot from Table 4 - File 



number 25. AR6051 (2213 UTC, 1/7/87). 
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number 26. AR6051 (2213 UTC, 1/7/87). 
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Figure 33. Scatter plot from Table 4 - File 

number 27. AR6054 (1508 UTC , 2/7/87). 
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Figure 34. 



Scatter plot from Table 4 - File 
number 28. AR6054 (1508 UTC, 2/7/87). 
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SHIP TRACI' SCATTER PLOT 




Figure 35. Scatter plot from Table 4 - File 

number 70. AR6113 (1526 UTC, 15/7/87). 
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Table 5 displays the differences in the ship track pixel 
means from the ambient pixel means. Table 5 also shows the 
conversion of channel 3 radiance to reflectivity. Table 5 
shows that there is a marked increase in values for channels 
1 and 3 in the ship track pixels compared to ambient pixels. 
Channel 1 ship track differences increased by 1.3% compared to 
ambient pixel means. Channel 3 ship track differences 
increased by 18.6% compared to the ambient pixels. One can 
also see that there is only a very slight increase in 
channel 4. This is again in support of the findings of 
Coakley et al. The percentage of increase in the values of 
reflectivity and radiance was found to be less in this study 
than that found in the results from Coakley et al., as seen in 
Table 6. Coakley et al., assumed that comparing the ratio of 
change in reflectivity for channel 1 to that for channel 3, 
removed the effect of fractional cloud cover which could 
influence the changes in the reflectivities noted above. The 
ratio of change in reflectivity for this analysis was 
0.18 +/-1.4 which is nearly half the value of 0.4 +/- 0.8 
found by Coakley et al. 
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TABLE 6 

COMPARISON TO COAKLEY et al. (1987) 


CHANGE IN 


% 


RATIO OF 


RADIANCE 


REFLECTIVITY 


CHANGE IN 


(mWm' 2 -sr' , -cm) 


INCREASE 


REFLECTIVITY 


CHANNEL 4 


CHANNEL 1 


CHANNEL 3 


CHN1 / CHN3 




COAKLEY et 


al. (1987) 




0.0+/-0.05 


1.6+/-0.7 


3.9+/-0.4 


0.4+/-0.8 




LUTZ 




<N 

rH 

• 

0 

1 

+ 

o 

• 

o 


0.34+/-2.6 


2.1+/-0.99 


0.16+/-2.6 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The primary focus of this thesis involved evaluating a 
large satellite data set to determine the number and 
geographical position of ship tracks that occurred during the 
period of study (1 - 19 July 1987). In total, 316 ship tracks 
were located and plotted. A statistical evaluation was then 
conducted on the radiative properties of the confirmed ship 
tracks. 

It was discovered that 22 of the 28 (76%) satellite passes 
evaluated in this thesis contained ship tracks. The 76% was 
shown to be significantly higher than the 5 to 10% expected in 
the earlier by Coakley et al. (1987) . It was shown that the 
AVHRR channel 3 was more favorable to the observance of ship 
tracks than channel 1. 

The large number of ship tracks found, provided a rich 
data set which was evaluated using the modified Coakley 
algorithm for automatic ship track detection. The statistical 
study concurred with the conclusions of Coakley et al., that 
ship exhaust increases the reflective properties of clouds. 
Table 4 indicated there is substantial variability of the 
cloud radiative properties, which show up clearly when 
examining the ship track and ambient pixels for 
channels 1, 3 and 4. When considering the means, all three 
channels showed an increase in the ship pixels compared to the 
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ambient pixels. The increase for channels 1 and 3 was 
significant (0.35 +/- 2.6% reflectivity for channel 1 and 

0.11 +/- 0.05 mWm‘ 2 sr' 1 cm for channel 3) and very slight for 
channel 4 (0.01 +/- 0.12 mWm' 2 sr' , cm) . Although all three 

channels showed increases in radiative properties, as in 
Coakley et al., the increases were less for this study. 

Many areas of research have yet to be explored in the 
study of ship tracks. Some of the areas for future work are: 

1. Investigation of weather patterns within areas of known 
ship track formation should be conducted to determine 
more precisely what factors encourage the formation of 
ship tracks. 

2. An effort should be made to discover the cause of the 
variability found in the reflectivity and radiance 
readings associated with the ship track. 

3. Work should be continued to improve the ability of the 
Coakley algorithm to process satellite data and auto- 
matically locate and analyze ship tracks. Additionally 
the algorithm should be developed to automatically create 
a plot of ship tracks found. 

4. Studies are still needed to determine if specific ship 
types can be associated with variations in ship track 
properties . 

The study of ship track formation will be a continuing 
source of exciting research in the years to come. Civilian 
and military use of the ability to monitor the movements of 
ships by evaluating their ship tracks is very promising. With 
the growing worldwide concern about our environment, the 
efforts expended on ship track formation processes will pay 
great dividends in our efforts to protect our fragile planet. 
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